home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / source / src / plane / _Point1.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-11-16  |  3.3 KB  |  164 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  _Point1.c
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16. #include <math.h>
  17. #include <ctype.h>
  18. #ifdef AMIGA
  19. #    include <LEDA/_Point1.h>
  20. #else
  21. #    include <LEDA/Point1.h>
  22. #endif
  23.  
  24. //------------------------------------------------------------------------------
  25. // Points 
  26. //------------------------------------------------------------------------------
  27.  
  28.  
  29. Point_rep::Point_rep()  { count=1; x = y = 0; w = 1;}
  30.  
  31. Point_rep::Point_rep(Int a, Int b) 
  32. { x = a; 
  33.   y = b;
  34.   w = 1; 
  35.   count = 1; 
  36. }
  37.  
  38.  
  39. Point_rep::Point_rep(Int a, Int b, Int c) 
  40. { x = a; 
  41.   y = b;
  42.   w = c; 
  43.   count = 1; 
  44. }
  45.  
  46. Point::Point()                  { PTR = new Point_rep; }
  47. Point::Point(Int x, Int y){ PTR = new Point_rep(x,y); }
  48. Point::Point(Int x, Int y, Int w){ PTR = new Point_rep(x,y,w); }
  49.  
  50. //Point::Point(vector v)          { PTR = new Point_rep(v[0], v[1]); }
  51.  
  52. /*
  53. double Point::angle(const Point& q, const Point& r) const
  54. {
  55.   double cosfi,fi,norm;
  56.   
  57.   double dx  = q.ptr()->x - ptr()->x; 
  58.   double dy  = q.ptr()->y - ptr()->y; 
  59.  
  60.   double dxs = r.ptr()->x - q.ptr()->x; 
  61.   double dys = r.ptr()->y - q.ptr()->y; 
  62.   
  63.   cosfi=dx*dxs+dy*dys;
  64.   
  65.   norm=(dx*dx+dy*dy)*(dxs*dxs+dys*dys);
  66.  
  67.   cosfi /= sqrt( norm );
  68.  
  69.   if (cosfi >=  1.0 ) return 0;
  70.   if (cosfi <= -1.0 ) return M_PI;
  71.   
  72.   fi=acos(cosfi);
  73.  
  74.   if (dx*dys-dy*dxs>0) return fi;
  75.  
  76.   return -fi;
  77. }
  78.   
  79.  
  80. Point Point::rotate(const Point& origin, double alpha) const
  81. { if (origin == *this) return *this;
  82.   segment s(origin,*this);
  83.   return s.rotate(alpha).end();
  84. }
  85.  
  86. Point Point::rotate(double alpha) const
  87. { return rotate(Point(0,0),alpha); }
  88.  
  89. Point Point::translate(double alpha, double d) const
  90. { double dx = cos(alpha) * d;
  91.   double dy = sin(alpha) * d;
  92.   return Point(ptr()->x+dx,ptr()->y+dy);
  93.  }
  94.  
  95. Point Point::translate(const vector& v) const
  96. { return Point(ptr()->x+v[0],ptr()->y+v[1]);
  97.  }
  98.  
  99. double Point::distance(const Point& p)  const
  100. { return hypot(p.ptr()->x - ptr()->x, p.ptr()->y - ptr()->y); } 
  101.  
  102. double Point::distance() const
  103. { return distance(Point(0,0)); }
  104. */
  105.  
  106.  
  107. int Point::operator==(const Point& p) const 
  108. { return ( (((ptr()->x) * (p.ptr()->w) - (p.ptr()->x) * (ptr()->w)) == 0) && 
  109. (((ptr()->y) * (p.ptr()->w) - (p.ptr()->y) * (ptr()->w)) == 0)); }
  110.  
  111.  
  112. ostream& operator<<(ostream& out, const Point& p)
  113. { out << "(" << p.X() << "," << p.Y() << "," << p.W() << ")";
  114.   return out;
  115.  } 
  116.  
  117. istream& operator>>(istream& in, Point& p) 
  118. { // syntax: {(} x {,} y {,} w {)}   
  119.  
  120.   int x,y,w;
  121.   char c;
  122.  
  123.   do in.get(c); while (in && isspace(c));
  124.  
  125.   if (!in) return in;
  126.  
  127.   if (c != '(') in.putback(c);
  128.  
  129.   in >> x;
  130.  
  131.   do in.get(c); while (isspace(c));
  132.   if (c != ',') in.putback(c);
  133.  
  134.   in >> y; 
  135.  
  136.   do in.get(c); while (isspace(c));
  137.   if (c != ',') in.putback(c);
  138.  
  139.   in >> w; 
  140.  
  141.   do in.get(c); while (c == ' ');
  142.   if (c != ')') in.putback(c);
  143.  
  144.   p = Point(x,y,w ); 
  145.   return in; 
  146.  
  147.  } 
  148.  
  149. /*
  150. int compare(const Point& a, const Point& b)
  151. { Int x = a.X() * b.W();
  152.   Int y = b.X() * a.W();
  153.   int s = sign(a.W()) * sign(b.W());
  154.   if (x == y) 
  155.   { x = a.Y() * b.W();
  156.     y = b.Y() * a.W();
  157.     if (x == y) return 0;
  158.    }
  159.   return (x < y) ? -s : s;
  160. }
  161. */
  162.  
  163.  
  164.